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> 1 Introduction 

Quantum systems in finite dimensional Hilbert spaces proved to be a good laboratory for inves- 
tigations of connections between quantum and classical characteristics of chaotic and integrable 
f — | dynamics. In particular statistical properties of quantum spectra are well established signatures 
of integrability on the classical level [HaaOOj. 

Recently we investigated a general procedure of attaining the classical limit for quantum 
systems with Hamiltonians defined as polynomials in generators of finite-dimensional compact 



Lie algebras |SK06j . We generalized to an arbitrary algebra g the procedure proposed in [GK98] 
and [GHKOOj , where a particular example of the SU3 algebra was used to illustrate connection 
between spectral statistics and properties of the limiting classical Hamiltonian system. In the 
above cited papers we discuss at length the physical relevance of model Hamiltonians defined 
in terms of generators of an arbitrary g, as well as the meaning of (various) classical limits of 
them, so here we only mention that they describe e.g., collections of many (say N) identical 
multilevel (or spin-like) systems. The classical limit is attained when N goes to infinity and 
is performed by increasing to infinity the dimension of the irreducible representation in the 
space of which the Hamiltonian of the system acts. In cases when the rank of the Lie algebra 
is greater than one, no unique way of attaining the classical limit via a sequence of irreducible 
representation with increasing dimensionality can be, a priori, singled out. This fact has some 
interesting consequences discussed thoroughly in the above cited papers (like e.g., different 
limiting statistics of the eigenvalues of the same system, but for different paths to the classical 
limit). 

In the present paper we want to study more closely spectral statistics of Hamiltonians 
defined in terms of elements of SU3 along sequences of irreducible representations mimicking 
the transition to the classical limit. Formally, the classical limit, is often treated as 'going with 
the value of the Planck constant to zero'. Such a statement taken literally is meaningless - we 
should specify with respect to value of which physical quantity the Planck constant becomes 
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negligible, so that the classical description of the system is justified. In the above mentioned 
situation of N M-level systems and a 'nonlinear' Hamiltonian in the form of a polynomial in the 
generators of su n the classical limit was obtained for N — > oo, and the 'effective Planck constant' 
(the quantity which approached zero in the classical limit) was scaled appropriately with the 
dimensions of representations. Such a scaling was motivated by purely physical considerations. 
In the present paper we show that for an arbitrary g there is a unique proper scaling relevant 
for the classical limit. This is one of the main result of the paper justifying our approach to 
study the outlined class of models in order to understand the transition to classical limit in 
well controllable models. 

In the following we will concentrate on Hamiltonians taking the form of second order poly- 
nomial in elements of SU3, inspired by the so called Lipkin Hamiltonian of a three-shell nuclear 
model, often invoked in various investigation of quantum- chaotic phenomena [LS90a, LS90b, 
WIC98J. We will be able to thoroughly investigate limiting spectral properties of the linear part 
of the Hamiltonian in all possible cases of attaining the classical limit and give some remarks on 
the properties of the full nonlinear Hamiltonian, which will be an object of our further studies 
in a forthcoming paper. 

Let us briefly sketch the contents of the remaining sections. In the following one we give the 
definition of the nearest neighbour distribution of eigenvalues of a hermitian matrix - the main 
object characterizing a quantum system investigated in the paper. Section [3] is concerned with 
rays of irreducible representations and an abstract notion of hermitian operator. In Section H] a 
limit theorem for simple operators, i.e., Lie algebra elements is given. Under some assumptions 
on the starting representation these operators will have Dirac statistics as limiting nearest 
neighbour distribution. A proof of this is given by estimating the number of different eigenvalues 
of the represented operator along the ray through the starting representation. As a corollary 
we obtain Dirac statistics for certain monomials in U(q). 

If we consider operators which are not given as monomials there is need for a suitable choice 
of rescaling. The already announced main result of Section is that under all possible ways of 
rescaling only the one used in [GHKOOJ and [SK06J gives interesting and physically meaningful 
results. For this we give physical arguments as well as mathematical facts. 

In the last section we deal with an explicit 'nonlinear' example. We consider the full Lipkin- 
Hamiltonian. By giving explicit formulas for the action of operators in arbitrary representation 
we are able to complete our statements from Section H] concerning the nearest neighbour statis- 
tics for the linear part for arbitrary choice of the sequence of representations and make some 
remarks about the structure of the matrix of the full Lipkin Hamiltonian in arbitrary represen- 
tation. 

2 Nearest neighbour distributions 

The statistics we deal with are the nearest neighbour distributions which are defined as follows 

Definition 1. Let A be an N x N hermitian matrix with eigenvalues, counted with multiplicity, 
X\ < • ■ ■ < Xn ■ 
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The nearest neighbour distribution of A is the Borel measure on R given by 
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ifx x = . . . = X N . 

In order to have a notion of limit of these distribution we make use of the following distance 
between measures which is common in statistics 

Definition 2. Let a, v be Borel measures on R>o of finite mass. The Kolmogorov-Smirnov 

distance of a and v is given by 



We will use the notion of convergence induced by this distance exclusively in text. Since 
weak convergence of measure is implied by the pointwise convergence of distribution functions 
the convergence in the sense of Kolmogorov-Smirnov also implies weak convergence. 



We fix a semi-simple, compact Lie group K and its complefixication G whose Lie algebras are 
denoted by t and g. We denote the universal enveloping algebra of g by U(g). 

By convention each considered representation of K will be irreducible and unitary and there- 
fore finite-dimensional. To have a simple notation we use the same letter for the representations 
of K and G and the induced Lie algebra representation. Moreover, this letter will be used for 
the induced representation on U(g). If p : K — > U(V) is such an irreducible, unitary repre- 
sentation, p : U(q) — ► End(V) is thus the induced representation on the universal enveloping 
algebra from the original p. 

There are numerous possibilities of sequences of irreducible representations of K, but using 
the results from [SK06J a physical meaningful limit is provided by rays of irreducible represen- 
tations. 

Definition 3. Let p : K — > U{V) be an irreducible representation of highest weight A. 

A ray through p is a sequence of irreducible representations (pk '■ U{q) — > End(Vk))km > 
such that pk has the highest weight k ■ A. 

We would like to emphasize that the irreducible representations of highest weight k ■ A are 
exactly the leading irreducible representations in the /c-fold tensor product of the representation 

P- 

Since we would like to treat the operators to be discussed as Hamiltonians of some physical 
systems, we need a representation independent notion of hermitian operators. 




(3) 



3 Setting 
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Lemma 4. There exists a unique M.-linear map f : U(g) — > U(q) satisfying 
• £ = forget 

. (<£ . . = cg-.-d MCj G W(j),cG C 
and /or a// £ G W(g) 

C f = £ (4) 
is hermitian, for each irreducible, unitary representation p : K End(V). 

Proof. First, one defines * on the full tensor algebra T(g) by the above requirements and checks 
by a direct calculation that is passes to the quotient. 

The statement about the image follows by another calculation from the defining prop- 
erties and the fact that p :U(g) — > End(V) is an algebra homomorphism. □ 

Thus, we may propose the following definition. 

Definition 5. An operator!; G U(q) is called abstract hermitian operator, if £ = £L The real 
subspace of abstract hermitian operators in U(q) is denoted by 7i. 

Remark 6. Note, that £ may not be an abstract hermitian operator although is hermitian, 
e.g., every £ in the kernel of p is represented as hermitian matrix. 



4 Simple Operators 

In this section the spectral statistics of simple operators, i.e., hermitian generators of a semi- 
simple compact Lie group, are considered along sequences of irreducible representations. The 
main result of this section is the following 

Theorem 7. If G is simple, rank G > 2 and p is an irreducible representation of highest weight 
\, such that A is in the interior of the Weyl chamber, then for all £ G it 

lim p Pk{0 = 5 (5) 

The strategy for the proof is as follows. First, £ is contained in the complexified Lie algebra 
t c of some fixed maximal torus T. We compare the dimension of the representation space 
with the number of different weights of the representation, since the weights evaluated at £ are 
the eigenvalues of Pk{0- The quotient of these numbers is shown to converge to zero, hence in 
the limit we will recover the Dirac measure for the nearest neighbor distribution of Pk{0- 

We start by giving an estimation the dimension of an irreducible representation in terms of 
the highest weight. 

Lemma 8. Let p be an irreducible representation with the highest weight A and A = Ylj \fj 
the decomposition of A into the fundamental weights fj . 
Then the dimension of p\ is bounded from below by 

n M. («> 

aen+,(A,a}>o x ' ' 
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where IT + denotes the set of positive roots and 5 — | Y] 



2 L^a<^X\~ 



a. 



Moreover, the number n\ of possible weights of p\ is bounded from above as follows 

n x < ordW- J](A, + 1). (7) 



Proof. The Weyl's dimension formula [W G99] reads 

(5 + A, a) _ -i-r f A i (A, a) 
+ (5, a) U 



Now, (A, a) > and (5, a) > for all positive roots a. Thus, the first inequality is clear. 

Starting from A we get all other weights by subtracting multiples of the roots. The lattice 
of roots is a sublattice of the lattice of weights, so we can reach every weight by subtracting 
multiples of the fundamental weights fj. 

There are at most rX/(-\? + 1) of such possible substractions which give positive weights 
and every weight is in the jy-orbit of a positive weight, which has at most \W\ elements. This 
proves the second inequality. □ 

Now, we give a proof of Theorem [71 

Proof of Theorem Let Q be the set 

Q := {a E Il + : (A, a) > 0} (9) 

and q := cardQ- We claim that Q = H+. Indeed, since A is in the interior of the Weyl chamber, 
we have 

(A,ay)>0 (10) 

for all simple roots <x,. But the positive roots are just positive linear combinations of simple 
roots, thus Q = U + . 

There exist exactly r := rankG simple roots and q > r because G is assumed to be simple 
and r > 2. If q = r were true, the root system could be decomposed into 1-dimensional pieces, 
which would contradict the assumptions. 

Combining the above two estimates, (jSJ) and (J7J), we obtain 

#differ. eigenvalues of p fc (Q ord(WQ II^Aj + 1) . Q .. r _ fl (u] 

^ K HaeQ {6,a) 

since the number of different eigenvalues is less then the number of weights of the representation. 
The right-hand side of (Hip converges to zero if k goes to oo which finishes the proof. □ 

Remark 9. Theorem UJ is still valid for G semi-simple provided that q > rankG. 
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5 Rescaling 



If we consider operators £ G H which are not 'homogeneous', e.g., £ = rji + r] 2 r] 3 , there is a 
question of the right scaling of the operators when the dimensionality of irreducible representa- 
tions increases. It will in general happen that the distance between the maximal and minimal 
eigenvalues of 772773 along a ray of representations grow much faster than the distance for 771. 
Thus, as will be shown, the lower degree terms vanish in the limit. It is therefore necessary to 
rescale the operators in dependence of the parameter k. 

Definition 10. Fix a basis £1, • • • ,£ n of g and identify U{g) with C[£i, •■,£«,]■ Let p : K — > U(V) 
be an irreducible representation and (pk)k<m be the ray through p. 

A rescaling map is a family of maps {r k :U(g) — > W(s))fceN given as 

r k (Q = (12) 

on the generators and extended multiplicatively, where each Sk € N* and linifc^oo s k = 00. 

In other words we put the elements in U(g) in normal ordering with respect to the chosen 
basis and substitute for each £j. Thus, rescaling is induced by the linear automorphism of 
g given by multiplication with s^ 1 . Although this automorphism is independent of the chosen 
basis, the resulting extension to a vector space automorphism of U(g) is not, since it depends 
on the identification of U{g) with C[Xi, . . . , X n ]. More precisely, it depends on the ordering of 
basis vectors. Two different scaling procedures can be termed 'natural' 

1. Rescaling by the dimension: := dimp^ 

2. Rescaling by the parameter: Sk '■= k for k > and So := 1. 

Our main goal in this section is to show that only the latter will give physically interesting 
results. 



5.1 A Physical Argument 

As indicated in the introduction and explained in details in [GK98[ IGHKOOl ISK06] , we want to 
treat the limiting procedure described in the previous paragraphs as corresponding to taking 
the classical limit of the physical system in question. Using the momentum map we associate 
to an irreducible representation p : K — > U(V) with highest weight A the coadjoint orbit K.X in 
t*. The classical limit is the performed exactly by increasing k to infinity along the ray through 
p the corresponding sequence of coadjoint orbits. In this way k — > 00 should correspond to 

0. 

We will show that k can be identified with Yc x by invoking a semiclassical argument to 
express the volume the phase space in terms of the Planck constans. We will show that the 
rescaling by inverse parameter is the only natural choice, since a first order operator should be 
measured in the natural scale, i.e., in terms of h, and their rescaling should be 

r k {Q = K 3 = ^{Q. (13) 

Operators of higher degree should then be rescaled accordingly, which implies rescaling by the 
parameter. To this end let us formulate 
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Proposition 11. Let p be an irreducible representation with highest weight X in the interior of 
the Weyl chamber. Then 

k = 1. (14) 

Proof. For the proof of (1141) recall that the volume of the phase space should be proportional 
to 

w (15) 

where / denotes the number of degrees of freedom in the classical phase space, i.e. the dimension 
of the coadjoint orbit through the highest weight. By a standard fact of the theory of compact 
Lie groups [] this number is equal to the number of positive roots, i.e., 

/ = # n + . (16) 

But the volume of the phase space K.X can be expressed in terms of the dimension of an 
irreducible representation. For this let 5 denote half of the sum of positive roots and set 
A' = A — 5. It is proved in [Kir68] that 

vol(AT.A) = dimp v . (17) 

Using Weyl's dimension formula it follows that 

= n = n M- («) 

, J - J - l (o, a) 

aen+ N ' ' aen+ x ' ' 



vol 



Therefore the volume of K.(kX) is proportional to k #n+ = k? . Comparing with the propor- 
tionality factor in ( fl5l) implies 

* = S (19) 

and completes the proof. □ 



5.2 A Mathematical Argument 

Although the reader may be already convinced by the above reasoning, we would like to give a 
more mathematical argument, which applies also to representations on the border of the Weyl 
chamber. 

We state the following lemma: 

Lemma 12. Let = ^ / a/S / G U{q) be given with £} H = and consider the ray (p& : K — > 

U(Vfc))fceN* through an irreducible representation p : K —>■ U(V) of highest weight X. 
Then 

\\p k ^H)\\End(V k ) <X)l^l C ( A ) |/LA;|11 ( 2 °) 

where the Cj(X) is a constant depending only on X and \\ ■ ||End(v fc ) denotes the operator norm on 
End(V k ). 
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Proof. We use the explicit construction of irreducible representations by Borel-Weil. For this 
let r := rankG and 

Sj = (s?,.. .,s%), j = l,...,r (21) 

denote a basis of the j-th fundamental representation. These are holomorphic sections in a 
holomorphic line bundle 

Lj -> G/B (22) 

where B is a Borel subgroup of G and L — G x x . C, such that %j : B —>■ C is the exponentiated 
character of the fundamental weight Xj. The irreducible representation with highest weight A 
is then given by the action on sections of the line bundle 

L = Lf Xl ® . . . ® L® Aj -> G/B. (23) 

By the theorem of Borel-Weil the tensors of the form 

S?®...®S£ r , (24) 

with Ii,...,I r multiindices of degree = Xj constitute a generating system of the space of 
sections. 

Without loss of generality we may assume that £i is represented by a diagonal hermitian 
matrix in every fundamental representation, whose spectral norms a bounded by a constant 
C£j . Since the operator norm is equal to the spectral norm, we wish to give an estimate for the 
maximal absolute value of an eigenvalue of £i in p\. 

But on the generating system of vectors given by ( 1241) the action is on each factor separately, 
so we have 

||p(6)ll <c 6 (A 1 + ... + A r )=:c 6 |A|. (25) 

For this recall, that £j acts as differential operator of degree 1, therefore we get the A/s as 
scalars and not in the exponent. 

Clearly, the same argument can be carried out for £ 2 > • • ■ >6n- Thus, we have the following 
estimate 

\\p k (Q\\<ck(Xi + ... + X r ) = ck\X\ (26) 

for all j — 1, . . . , n with the constant c := maxj eg.. 
Now, consider 7 = J2i ^iX 1 . Then 

||Pfc(7)l|End(V fe ) < ^2 l a/ ill^ fc (£OIILd(Vy ' • • • ' ll^(^n)llEnd(F fc )- 

(27) 

1 

Using the estimates given by (126]) and Lemma [H we see that 

lk(7)lkd ( v fe )<El^i A;|/L ( c i A l) |/ ' ( 28 ) 

1 

which completes the proof. □ 

Remark 13. If £ G g, £' = ^ and pi(£) 7^ 0, then the asymptotic growth of the norm is given by 

\\ Pl (0\\~c-k, (29) 

where c is a constant, which depends on A and ||pi(£)ll- This can be seen by a short calculation 
for the action of £ on the tensors in (pTj) . 
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Recall that for a highest weight A the set Q is defined as Q = {a G H+ : (a, A)} and 
q = #Q. We use the Lemma and the remark to prove the following theorem 

Theorem 14. Consider the ray (p m : K — > U(K„)) me N* through an irreducible representation 
p : K — > \J(V) of highest weight A and assume q > 2. 

Lei £h = £ + Xl|/|>i a /" / £ W(g) fre given, such that £' H = an d £ G wift ||pi(C)ll 7^ 0- 
T/ien i/ie nearest neighbor distributions of rescaled a/ong' t/ie ra?/ agree wit/i t/ie /imit o/ t/ie 
nearest neighbor distribution of if the res caling factor sk grows at least like k 1+e . 

Proof. We claim that ||pfc(^fc(X]|/|>i) a i'" / )llEnd(y fc ) - ► 0, as k — > oo. Making use of (1251) and the 
growth of Sfc we obtain 



End(V fe ) S 



k^clXl^ k^ , , 

* —Fir- * c fcOT' ( 3 °) 



where c is a constant. On the other hand by the remark above Pk(j~C) has norm bound from 



below by some constant. Thus, \\pk(rk(0) ~ p( r fc(6/))||End(v fc ) -> as k -> oo. □ 

Analogously, one would like to treat the case of "underscaling" , i.e., if grows slower than 
k l ~ e . In that case we would expect the limit of the nearest neighbour distribution to depend 
only on the top homogeneous part of While this is the case in all known examples, we can 
give no proof here. The problem being that a lower bound for the growth of the norm in the top 
homogeneous part is needed. If the norms of the top homogeneous part behave asymptotically 
like in ( 1281) then with an analogous argument like in the proof above one can show convergence 
to the nearest neighbour distribution of the top homogeneous part. 



5.3 A remark on the commutativity 

Let us briefly discuss rescaling from the semiclassical viewpoint, by pointing that the postulated 
rescaling ensures the commutativity of the dynamical variables in the semiclassical limit. For 
simple operators let's consider rfc(£i£ 2 ) = = t-^2^i + t-^i,^]- As the commutator is a 

S k S k S k 

linear combination of £j, the respective term will vanish as k goes to infinity. A similar reasoning 
can be applied to arbitrary homogenous polynomials in ^ . So, in this narrow sense the amount 
of commutativity is increased, although the operators themselves are non-commutative. We 
like to think of this as a sign of the semiclassical nature of the limit. 



6 The Lipkin-Hamiltonian 

In this section we will make some remarks about more complicated operators on SU3, i.e. simple 
"polynomials" in the universal enveloping algebra of $1(3, C). 

Recall that the fundamental irreducible representations pk of the group SU n are given as 
the natural representations of SU n on the vector spaces 

k 

/\C n for Jfe = l,...,n-1. (31) 

i=l 
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By the theorem of Peter and Weyl every irreducible representation is realized as a subrepre- 
sentation of the natural SU n representation on L2(SU n ). It is known, what the generating 
functions for them are: let g G SU n be the matrix 

'9u ■■■ sO 

i ; , (32) 

\g n l ■ ■ ■ Qnn J 

then the vector space of the irreducible representation with highest weight A = Y^k=i ^kfk, 
where fk denotes the highest weight of pk is given as the linear span of the fc-homogeneous 
polynomials in the minors of g of the appropriate degree taken from the first columns only. 
This means, we take the homogeneous polynomials of degree Ai in the 1 x 1-minors of the first 
column, the homogeneous polynomials of degree A2 in the 2 x 2-minors of the first two columns 
and so forth, and then take every possible multiplicative combination of them. 
For SU 3 we have to consider only the 1 x 1-minors 

xi := g n , x 2 := #12, x 3 := 5-13 (33) 

and the 2 x 2-minors 

Vl ■= 921932 - 931922, V2 ■= 93l9l2 ~ 9ll932, V3 ■= 911922 ~ 92l9l2- (34) 

A generating set for the representation subspace of L 2 (SU 3 ) corresponding to the highest weight 
A = (Ai, A2) is then generated by the polynomials in the x« and y^ which are homogeneous of 
degree Ai in the Xi and homogeneous of degree A 2 in the yj. Unfortunately, these are not always 
linearly independent due to the relation 

xm + x 2 y 2 + x 3 y 3 = 0. (35) 

It can be shown, that this is the only relation among the variables, i.e., the irreducible 
representations are given as the bi-graded subrepresentations of the representation on the C- 
algebra 

R:= C[x l ,x 2 ,x 3 ,y 1 ,y 2 ,y 3 ]/(x 1 y 1 + x 2 y 2 + x 3 y 3 ). (36) 

The necessary details for calculating these relations for arbitary SU n can be found in [MS05J 
and [DRS74j . In the first book there is actually a nice interpretation in terms of Grobner bases. 

It was shown in the last section that elements from su.3 will have spectral statistics converging 
to the Dirac measure in the interior of the Weyl chamber. We want to study a little more 
complicated object, called the Lipkin-Hamiltonian in nuclear physics (cf.[GHK00j). It is an 
abstractly hermitian operator in T(s [3(C)) given by 

pipkin = aT 3 + b^2 S ij> ( 37 ) 

where T 3 = diag(l, 0, —1) and S^- are the matrices with 1 at position and elsewhere and 
a, b denote some real constants. For b = 0, we already know that the spectral statitstics in the 
interior of the Weyl chamber converge to zero. Before we go any further, we will first calculate, 
how these operators act. 
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We start with the example of 5*12. It acts on a vector v by j| | exp(Si 2 t).v . Since we are 
acting on functions, the action on the argument is given by multiplication with the inverse. 

/ #11 #12 9i3\ / 9n ~ tg 2 i #12 ~ tg 2 2 #13 ~ tg 23 \ 
exp(Si 2 ty l g 21 g 22 g 23 \ = \ g 2l g 22 g 23 . (38) 

\#31 932 933/ \ 931 g32 #33 / 

Thus, we get the following substitution rules for the arguments Xi and yc 

Xi = gn ' * #11 - tg 21 = xi - tx 2 
x 2 = g 2 i i * x 2 
%3 = 9si i-> 2:3 

Vl = 921932 - 93\g22 >-> 2/1 

2/2 = Stasia - #n#32 i-> #31 (#12 - ^12) - (#11 - tg 2 i)g32 = 92 + tyi 

93 = #21#32 - #31#22 l~> (#11 ~ ^#2l)#22 - #21 (#12 ~ ^#22) = #3 



(39) 



Applying this to the basis elements x J ll x 2 2 x 3 i3 y 1 1 y 2 2 y 3 3 the action is given by 

exp{S X2 t).x?x%xfy\^y b i = -o^-^+^yf yfyj 

Thus, the S'12 action is given by the action of the differential operator 



d 
dt 



t=o 



(40) 



9 d 

-x 2 — + yi — 41 
dxi d# 2 

on the bi-homogeneous polynomials. In fact, analogously we can calculate that the action of 
Sij is given by 

- x >l, +v W/ (42) 

The actions of T 3 corresponds to the differential operator 

d d d d , AnS 

~ x ±7> hx 3^ \-9it, 2/3^—- (43) 

dxx dx 3 dyi dy 3 

Remark 15. The reader should be aware of the fact that in the chosen basis operators which 
are abstract hermitian, e.g., S'12 + S21, are not represented as hermitian matrices. 

We begin our discussion by showing that T 3 has the Dirac measure as limiting spectral 
statistic in any case. Recall, that we proved this only for the interior of the Weyl chamber in 
Section [4J 

Lemma 16. Let A be the highest weight of a non-trivial irreducible representation of SU 3 and 
assume that A lies on the border of the Weyl chamber of SU 3 , i.e. A = (Ai,0) or A = (0, A2) 
and A ^ (0,0). 

Then the spectral statistics of the operator T = diag(l, 0,-1) G sls(C) along the irreducible 
representations p*^ m \ converge to the Dirac measure for m — > 00. 



11 



Proof. For this we consider only the case A 2 = 0, the other case being analogous. 
The element T operates on a polynomial x^ 1 x a 2 2 a^ 1- " 1-02 by 

T.x a l 1 x a 2 2 xl 1 - ai ~ a2 = {- ai + X 1 -a 1 - a 2 )x a l 1 x a 2 2 x^ ai ~ a ' 2 . 



(44) 



So the eigenvalues of T are given by —2a! — a 2 + Ai with < a±, a 2 and a\ + a 2 < Ai. We see 
for example that the zero eigenvalue has increasing multiplicity if Ai increases, because 



= 2 • + Ai - Ai = 2 • 1 + (Ai - 2) - Ai = 2 • 2 + (Ai - 4) - Ai 



(45) 



Let us count the number of different eigenvalues now. This number is the number the fibers 
of the map (ai,a 2 ) i— > 2ai + a 2 , which is 2Ai — 1. To see this, just note that (ai,a 2 ) is in the 
same fiber as [a\ + l,a 2 — 2). So every fiber can be represented by an element of the form 
(0, a 2 ) or (ai, Ai) and these are exactly 2Ai — 1 elements. 

But the dimension of the representation on the space on homogeneous polynomials of degree 
Ai is |Ai(Ai — 1). Thus the quotient of the number of eigenvalues and the dimension is 

2Ai - 1 4 2 

lA^Ax-l) " aT^T " A^Ai-1) (46) 

and converges to zero if Ai — > oo. Thus, we have proved the lemma. □ 

Our basic input for the Lipkin Hamiltonian is the explicit action of the Sfj. These act by 
sparse matrices as we are going to discuss now. 

Analogously, the action of other SV,- can be calculated. We give the action in a tabu- 
lar form here, where we just write [. ..] := [a 1? a 2 , a 3 , &i, b 2 , b 3 ] instead of the rather clumsy 

xTx a 2 2 xfy^y h 2 2 y h i. 



5*12 
5*13 
5*21 
S 2 3 
5*31 
5*32 



-ai[ai - 1, a 2 + 1, a 3 , &i, b 2 , b 3 ] + b 2 [ai, a 2 , a 3 , bi + 1, b 2 — 1, 6 3 ] 
-ai[ai - 1, a 2 , a 3 + 1, &i, & 2 , b 3 ] + b 3 [a x , a 2 , a 3 , &i + 1, 6 2 , o 3 - 1] 
-a 2 [ai + 1, a 2 — 1, a 3 , &i, & 2 , 6 3 ] + h[a 1: a 2 , a 3 , &i - 1, b 2 + 1, 6 3 ] 
-a 2 [ai, a 2 - 1, a 3 + 1, &i, b 2 , 6 3 ] + 6 3 [ai, a 2 , a 3 , 6i, 6 2 + 1, b 3 - 1] 
-a 3 [ai + l,a 2 ,a 3 - l,6i,6 2 ,6 3 ] + b^, a 2 , a 3 , b x - l,b 2 ,b 3 + 1] 
-a 3 [ai, a 2 + 1, a 3 - 1, 6i, 6 2 , & 3 ] + b 2 [a u a 2 , a 3 , 6i, 6 2 - 1, b 3 + 1] 



(47) 



Applying every SV,- twice we get 



r-2 
°12 



r-2 
°13 



r-2 
°21 



[. . .] h^ai(ai - 1) [ai - 2,a 2 + 2, a 3 , 6i, & 2 , & 3 ] 

+ 6 2 (6 2 - 1) [ai, a 2 , a 3 , 6i + 2, 6 2 - 2, 6 3 ] 

- 2ai6 2 [ai - 1, a 2 + 1, a 3 , 6i + 1, 6 2 - 1, 6 3 ] 
[. . . ] h^ai(ai - 1) [ai - 2, a 2 , a 3 + 2, 6i, & 2 , & 3 ] 

+ 6 3 (6 3 - 1) [ai, a 2 , a 3 , &i + 2, & 2 , &3 - 2] 

- 2a!0 3 [ai - 1, a 2 , a 3 + 1, 6i + 1, b 2 , b 3 - 1] 
[. . . ] ^a 2 (a 2 - 1) [ai + 2, a 2 - 2, a 3 , 6i, 6 2 , & 3 ] 

+ 61(61 - 1) [ai, a 2 , a 3 , 6 X - 2, 6 2 + 2, & 3 ] 

- 2a 2 0! [ai + 1, a 2 - 1, a 3 , &i - 1, 6 2 + 1, 6 3 ] 



(46 
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From this we can read of the matrix representation of the Sfj on the vector space of bi-graded 
polynomials of bidegree (Ai,A2) in R, where we choose the representatives [a±, 0,2, 0,3, 61, &2, 63] 
with 0363 = 0, i.e. we replace every X3JJ3 by —xiyi — X2i/2- We see immediately that Sf 2 and 
S21 do not increase a 3 or 63, so the matrix representation of Sf 2 and S 2 i has at most 3 entries 
not equal to zero in each column. The other SV, have two summands which may increase either 
<23 or 63, so they have at most 5 non-zero entries in each column. Thus, the sum Yli^j Sfj nas 
at most 2 • 3 + 4 • 5 = 26 non-zero entries in each column, regardless of the bi-degree (Ai, A2). 
Each of these entries has absolute value bounded by max{Ai, A2} 2 . 

So, the it possible to compute the Lipkin Hamiltonian numerically and get graphical repre- 
sentations as in the following figure: 



Full Lipkin Hamiltonian, d=728 



Lipkin 
Poisson 
GOE 



1.5 - 



0.5 



Width 



Here, the nearest neighbor statistics are drawn for the Lipkin Hamiltonian in the irreducible 
representation corresponding to Ai = A2 = 8, which is 728 dimensional. 
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